
Calhoun 

iniQiuiic^iul Ar{hiv« of tilt Mil vdl Poii^roduiit School 


Calhoun: The NPS Institutional Archive 
□Space Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


1994-09 

A numerical study of wind forcing effects on 
the California Current System 


Vann, James R. 

Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/43035 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 

Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


htt p://w ww. n ps. e du/l ib ra ry 


Callwuo is the Naval Postgraduate School's public access digital repository for 
research mate rials and institutiional publicatkios created by the NPS community. 
Calhoun is named for Professor of Mathematics Guy K. Caftiouo, NPS's first 
appointed — and published — schoteily author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 Univefsity Circle 
Monterey, California USA 93943 






AD-A286 205 


NPS-OC-94-006 

FOSIGRADU^ Sai(X)^ 
Monter^, California 



THESIS 


A NUMERICAL STUDY OF 
WIND FORCING EFFECTS ON THE 
CALIFORNIA CURRENT SYSTEM 

by 

James R. Vann 
September, 1994 

Thesis Advisor: Mary L. Batteen 


^>proved for public release; distribution is unlimited. 


Prepared for: 

Office of Naval Research 
800 N. Quincy Street 
Arlington, VA 22217-5000 


n - 




' 94-35200 









NAVAL POSTGRADUATE SCHOOL 

CA 93943 


Rear Admiiil Thomas A. Mercer 
SiqjerinteiideBt 


This thesis was p rep ai ed in conjunction with r esear c h sponsored in part by 
the Office of Naval Research, 800 N. Quincy Street, Arliitgton, VA 22217- 
5000. 


Rqnodoction of all or part of this report b airthorixed. 


Released by: 









REPORT DOCUMENTATION PAGE 


Pom A p proved (M4B No. 0704 


1 




lalamniaB OpMaoM aad Raaom. 1313 Jaflanaa Dava Hi^airv. Mw ilOa. AitogWa. VA 223024102, la *a Otltaa at ItaaaMMal aad Badri. 

PaaaiwoA KiiaiaiBa Biaiaat (OTOAOltt) Wadaiaiaa DC 3QKB. 

1. AGENCY USE ONLY deaw MMl) 2. R£K»T DATE 3. REPORT TYPE AND DATES COVERED 

Sapiaadiar 1994 Master’s UteSif 

4 TITLE AND SUBTITLE A NUMERICAL STUDY OF WIND FORCING 
EFFECTS ON THE CALIFORNIA CURRENT SYSTEM (U) 

S. FUNDING NUMBERS 

6. AUTHOR(S) Vtnn, James R. aa nn^rtina wito Mary BaMaa 

7 PERFORMING ORGANIZATION NAME(S) AND ADDRESS(BS) 

Naval PoMfraduate Sdiool 

Monterey CA 93943-5000 

8. PERFORMING 

ORGANIZATION 

REPORT NUMBER 
NPS-OC-94-006 

9 SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(BS) 

OAoe of Naval Reaaarch 
•OON. QiMBcy St. 

AfteWtoa. VA 22217-SOOO 

10. SPONSOBUNGHdONITORlNG 
AGENCY REPORT NUMBER 

11 SUPPLEMENTARY NOTES The views expccaaed in this diedia are thoae of die amhor and do not reflect I 

the offldal policy or pototion of die Department of Defense or die U.S. Govemmem. | 

1 12a. DISTRIBUTION/AVAILABILITY STATEMENT 

1 Approved for public release; distribution ia unlimited. 

12b. DISTRlBimON CODE I 

___ 


13. ABSTRACT 0m»mum 200 wor4$) 


A te^MwotMioa, te eJ. priwii^i'r inirtna doom aodii i* Mad to nutoiaB i M po—t of aa idaafard, flat-lMaom, 
aaaiani liwwitoq' ooaaaic rafiaa oa a bato-plaaa to ciiaMtalogical avarapa (19S0-19t9), iadividual yaar, aad aaiWpla yaar wiad 
loraai. Tka focaa of toia aaidy ia toa Caiifataia Canaat Syalaa aloag *a ooaalal ragioa, froto 3S* N to 47 Jf* N, off Ok Waat 
Coaai of No(* Aaaarica. 

Two typaa of axparaaaaai are onadactod. The firat qrpe forcea dw andel froaa real wdh rJimalnlngical, 1981, aad 1983 
ai n ad J y wiada to aaiMiaa dto gaaaraliaa phaaa of fcaairai auch aa cumata, upwaOiag, aieaadefa, addiaa, aad filamwaa. The 

ia toe frai type of oapariaMato, toa foBowiag faaliiwn are o haar v ad; a pokavard ooaatal auikoe curraot aear the itait aad 
aad of each yaar, aa aq aa torw a rd auikoe correal, a polearard uadercuneot, upweDiog, awaadera, aad addiaa. la die aecood type 
of inpariaiaala, aMaadare aad addiaa were already preaeal at toe alait of die eaqieriaieat. Ia additioo to die featurea obaerved 
doriag toe firat type of ea peri a iaBt, filaiaeata are ge aeta to d. The reauBa aiippoit die hypolheais diat wind forciiig ia aa importaiit 
atocbaaiaaa for toe geaeralioa of aaaay of dw obaerved featurea ia the Califoraia Curreat Syatem. 


14. SUBJECT TERMS Primitive equation model, California Current System, eddies, 
filamenii, jets, undercurrett, wind forcing, currents. 

IS. NUMBER OF 

PAGES 

235 

16. PRICE CODE 

17. SECURITY CLASSin- 

18. SECURITY CLASSIFI- 

19. SECURITY CLASSffl- 

20. LIMITATION OF 

CATION OF REPORT 

CATION OF THIS PAGE 

CATION OF ABSTRACT 

ABSTRACT 

Unclamfied 

Unclassified 

Unclassified 

UL 


NSN 7S4041-2SO-5500 Standard Form 298 (Rev. 2-89) 


Pnacribed by ANSI Sid. 239-18 


1 













































Approved for public release; distribution is unlimited. 


A NUMERICAL STUDY OF 
WIND FORCING EFFECTS ON THE 
CALIFORNIA CURRENT SYSTEM 


by 

James R. Vann 

Lieutenant, United States Navy 
B.S., United States Naval Academy, 1984 

Submitted in partial fulfillment 
of the requirements for the degree of 

MASTER OF SCIENCE IN 

METEOROLOGY AND PHYSICAL OCEANOGRAPHY 

from the 


Author; 


^}proved by: 


NAVAL POSTGRADUATE SCHOOL 
September 1994 



Robert H. Bourke, Chairman, 
Department of Oceanognqrfiy 


m 













iv 





ABSTRACT 


A high-resolution, multi-level, primitive equation ocean model is used to 
examine the response of an idealized, flat-bottom, eastern boimdary oceanic 
regime on a beta-plane to climatological average (1980-1989), individual year, 
and multiple year wind forcing. The focus of this study is the California Current 
System along the coastal region, from 35* N to 47.5* N, off the West Coast of 
North America. 

Two types of experiments are conducted. The first type forces the model 
from rest with climatological, 1981, and 1983 monthly winds to examine the 
generation phase of features such as currents, upwelling, meanders, eddies, and 
filaments. The second type continues the forcing from the previous years to 
examine the maintenance of these features. 


In the first type of experiments, the following features are observed: a 
poleward coastal surface current near the start and end of each year, an 
equatorward surface current, a poleward undercurrent, upwelling, meanders, and 
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I. INTRODnCTIQN 


The California Current System (CCS) is aua eastern 
boundary current (EBC) off the west coast of the United 
States of America. A primitive equation numerical model 
developed by Batteen sa aJL. (1989) is used to simulate the 
currents, eddies and filaments in the CCS. In this study 
monthly winds from the Buropecm Centre for Medium-Range 
Weather Forecasts (ECMWF) global data set from 1980 to 1989 
(which include El Nifio events) are used to force the model. 

The thesis is orgamized as follows: Chapter II 
describes the region modeled, the El Nifio-Southem 
Oscillation (ENSO), and some previous modeling work in this 
region. The numerical model which was used here along with 
the experimental setup are presented in Chapter III. The 
results of the various model simulations are presented in 
Chapter IV. Chapter V summarizes the work, discusses naval 
relevance, and presents recommendations for future research. 
Comparisons of model simulations with availeJsle observations 
are shown in both Chapters IV and V. Six appendices follow 
which contain the detailed results of the barotropic cuid 
baroclinic instability amalyses for the main experiments. 
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II. BACXGRODIID 


A. BRIEF DESCRIPTION OF REGION 

The west coast of North America serves as an eastern 
boundary for the circulation of %fater in the North Pacific 
Ocean. The entire flow along the coast is referred to as 
the CCS, which is a classical subtropical EBC system. 

The CCS consists of three separate component flows, each 
of which varies both spatially and seasonally (see Figure 
2.1): an equatorward surface current, a poleward 
undercurrent and a poleward surface current adjacent to the 
coast (Tisch fit al., 1992). This latter current ceui be a 
surface manifestation of the undercurrent or a separate 
current. 

The predomincmt flow is an equatorward surface flow 
called the California Current (CC). The second basic flow 
is the poleward undercurrent which has been called both the 
California Undercurrent (CUC) (e.g., by Batteen fit al.. 

1989) and the California Counter Current (CCC) (e.g., by 
Hickey, 1979}. Thi^ is a relatively weak subsurface flow 
(around 5 cm s*^). 

The third basic flow is the Inshore Current (IC), which 
is known as the Davidson Current (DC) north of Pt. 
Conception. This is also a relatively weaJc poleward flow 
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found at the surface and near the coast. There are other 
ICs known as the Southern California Countercurrent (SCO to 
the south, and the Southern California Eddy (SCE) inshore of 
the Channel Islands within the California Bight (Hickey, 
1979), which are both equatorward of the model region (see 
below) to be studied. 

The CC is a broad, slow (~10 cm s'^) flow which can 
extend offshore around 1000 km (Hickey, 1979). Superin^osed 
on this mean flow are highly energetic mesoscale eddies and 
meandering jets which Ceui have peak velocities near 80 cm 
s~^ (Mooers and Robinson, 1984). Figure 2.2 depicts these 
features. 

A section of the CCS from 35® to 47.5® N emd extending 
12.5 degrees offshore from the coast will be the model study 
area and is depicted in Figure 2.3. Major observational 
progreuns which have sampled this region are quite numerous 
and include: California Cooperative Oceanic Fisheries 
Investigation (CalCOFI) , Coastal Ocecui Dynamics E3q}eriment 
(CODE) I and II, Large-Scale West Coast Shelf Experiment 
(nickncuned SuperCODE), Coastal Upwelling Experiment (CUE) I 
and II, Oceom Prediction Through Observation, Modeling and 
Analysis program (OPTC^), and Central California Coastal 
Circulation Study (CCCS). Figure 2.4 depicts some of the 
areas of observations. 
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B. CLZMATOLOQZCAL WINDS 


The North Pacific Subtropical High dominates the 
atmospheric circulation off the west coast of the United 
States of America, and varies in position and strength by 
season (Huyer, 1983). In the winter (Jcuiuary) the high Is 
about 1020 mb and centered near 30** N, 135** W (Figure 2.5a). 
In the summer (July) the pressure Increases to 1025 mb euid 
the center moves northward to 38® N, 150® W (Figure 2.5b). 
Over the continental west coast the typical pressure Is 1020 
mb In the winter, and 1015 mb during the summer, the latter 
In response to a thermal low (of 1005 mb) In the central 
valley of California (Huyer, 1983). This results In a 
tighter pressure gradient against the west coast during the 
summer. As this gradient tightens, stronger winds can 
result. 

Propagating atmospheric disturbances (Halllwell and 
Allen, 1987), along with other atmospheric mesoscale 
phenomena (Huyer, 1983), further conpllcate the wind pattern 
on a daily, as well as on a seasonal basis along the coast. 
The winds are also affected by coastal atmospheric boundary 
layer processes within 100 to 200 km of the coast (Halllwell 
2 Uid Allen, 1987) . 

The oceanic response to summer cmd winter wind patterns 
are most notable near the coast. The eguatorward winds of 
summer can produce nearshore upwelling, while the divergent 
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wind flow of winter cam force downwelllng in some areas 
(Rlenecker and Mooers, 1986). 

The shifting of the Pacific High, along with the 
Aleutian Low moving to the southeast, Ceui cause the winds to 
diverge near Cape Mendocino (located at edsout 40** N) during 
the winter. North of the Cape, northerly winds are 
encountered as the wind blows aroiind the eastern side of the 
Aleutian Low, while south of the Cape, southerly winds 
result as wind flows around the North Pacific High. 

The summer wind pattern is established by the 
interaction between the North Pacific High amd the thermal 
low in the southwest United States (Nelson, 1977; Halllwell 
and Allen, 1987). Winds are predominately eguatorward on 
both sides of Cape Mendocino during the summer (Nelson, 

1977). 

Climatological wind stress fields, shown in Figure 2.6 
for the months of January, April, July, amd October, depict 
regions of maximum wind stress in the CCS as shaded areas. 
The strongest wind stress is located off Point Conception in 
;^ril amd migrates northward until it is at a maocimum in 
July off Cape Mendocino. The average summer wind stress 
close to the shore is eguatorward amd is thus favorad>le for 
upwelling at the coast (Halllwell amd Allen, 1987; Nelson, 
1977; Wickham fit » 1987). 
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C. BL NifiO-SOUTHEBN OSCILLATION 

One large-scale two way coupling of the atmosphere and 
ocean has been studied for mainy years aind Is referred to as 
the El Nlfio-Southern Oscillation (ENSO) . El Nliio refers to 
the warming of Eastern Pacific oceem waters, while the 
Southern Oscillation is an atmospheric pressure fluctuation 
between the Indicui Ocean and the Eastern Tropical Pacific 
Ocean (Philcuider, 1990). Figure 2.7 illustrates key 
properties of the teleconnection including upper troposphere 
and surface pressure euiomalies. Note the shift in the 
Eastern Pacific Siibtropical High during cm ENSO event. 

Explanations for this phenomena involve both cm oceanic 
and atmospheric teleconnection. Rienecker and Mooers (1986) 
explain that the oceanic connection occurs when eastward 
propagating equatorial Kelvin waves turn poleward upon 
inpacting the west coast of the Americas. In contrast, the 
atmospheric teleconnection links the tropical sea surface 
tenperature (SST) anomalies to the Aleuticm Low by momentum 
transfer through am intensified Hadley Cell circulation 
(Bjerknes, 1966; Phllamder, 1990). A stronger poleward wind 
stress subsequently develops at the coast near 36** N when 
the Aleutiam Low deepens amd moves to the southwest. 
Simultaneously the North Pacific High weakens amd moves 
offshore. 

The El Niho of 1982 - 1983 v^as one of the strongest on 
record. The Aleutiam Low was much deeper than normal during 
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the winter and wind stress anomalies were at a maxixnum in 
the region between 35** N and 45** N (Trenberth sL al* # 1990). 
Anomalous conditions were observed from San Diego to 
Vancouver Island (Rienecker and Mooers, 1986). There were 
two anomalous peaks in the SST, the first occurring in 
February to March of 1983 off northern California, while the 
second one was evident in September to October of 1983 near 
Point Conception (Rienecker and Mooers, 1986), the latter 
equatorward of our model domain. 

D. NDKERZCAL MODEL STUDIES 

There have been niimerous models used over the past few 
decades to model EBCs. particularly the CCS. Extensive 
reviews Ceua be foxmd in Allen (1980), Chelton (1984), and 
O'Brien gt al* (1977). No attenpt is made to provide a 
conprehenslve review of the models here; instead, a sximmary 
of modeling work that has led directly to this study is 
presented. 

Early work Included that of Pedlosky (1974) and 
Philander cuid Yoon (1982) who used steady wind stress cuid 
transient wind forcing, respectively. Carton (1984) euid 
Carton and Phllemder (1984) investigated the response of 
reduced gravity models to realistic coastal winds. At 
periods exceeding 100 days the alongshore currents begcui to 
weaJcen and disperse as a series of alternating jets. 
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Another series of experiments were done by McCreary 
al . (1987) utilizing a linear model with both transient amd 
steady wind forcing with Laplaclan diffusion In which no 
eddies or fllaunents developed. This work was extended using 
a full primitive equation model by Tlelklng (1988) who 
obtained similar results. 

Batteen al* (1989) used a primitive equation, multi¬ 
level model with blharmonlc heat and momentum diffusion 
which had both steady forcing constauit throughout the domain 
and meridionally varying forcing. In both cases, when the 
barocllnlc shear between the surface current and 
undercurrent became strong enough, eddies and filaments 
developed. This was the first SBC model to simulate 
currents, eddies and jets In the CCS. 

The work of Batteen (1989) was extended by 

Mitchell (1993) cuid Haines (1994). Mitchell (1993) forced 
the model with different sets of averaged winds and obtained 
reasoncd)le results utilizing a set of seasonally varying 
full climatological wind fields. Haines (1994) enlarged the 
model domain, Inqproved resolution to 8 x 13 km euid used 
Intercuinual vice steady or seasonal winds. 

This effort builds on the model used by Batteen si. al . 
(1989) and the work of Mitchell (1993) and Haines (1994). 

The conputer code for the model has been modified to Include 
a barotroplc component, while the model resolution has been 
refined to 8 x ll km grid spacing. Seasonal auid Interemnual 
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wind forcing was used to force the model. Model simulation 
times have been expanded from one year, as in Mitchell 
(1993) and Haines (1994), to two to four years for each 
model run. The enhanced simulation times should allow 
analysis of the maintencuice of currents and eddies as well 
as the generation of these features. 
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Longitude (degrees west) 


Figure 2.1 Generalized circulation schematic of the 
California Current System (CCS). The broad, slow surface 
eguatorward California Current (CC) overlies the poleward 
California Undercurrent (CUC) along with the Inshore Current 
(IC), also known as the Davidson Current (DC), Southern 
California Countercurrent (SCC), and Southern California 
Eddy (SCE). 
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Figure 2.2 Section of California Current System with 
schematized positions of jets and synoptic-scale eddies; C, 
cyclonic eddy; A, anticyclonic eddy. Arrows Indicate 
principal surface flows; solid arrows are definite and open 
arrows are inferential. Also shown is a qualitative 
analysis of satellite Infrared image gray shades for sea- 
surface ten 5 )erature (SST) predominant patterns on l August 
1982: coolest SST, crosses; warmest SST, dots (from Mooers 
and Robinson, 1984). 


12 

























































































































































































































Latitude (degrees north) 


50 - 


47.5 


• • • • 


Columbia River 


• • 


• • 


42.5 


• • 


Cape Mendicino 


• • • • 


37.5 


• • • • 


Point Conception 

t • • • 


32.5 


• 


30 
-141 


I] 


Longitude (degrees west) 


Figure 2.3 Domain of model off the western coast of the 
United States of America. Domain bounded by 35® to 47.5® N 
and from the coast to 12.5 degrees offshore. Model covers 
an area 1024 km (1408 km) In the cross-shore (alongshore) 
direction and the cross-shore (alongshore) horizontal 
resolution is 8 Icm (11 Ion) . 








50'W 



Figure 2.4 Map of a portion of the west coast of North 
America, showing the location of some of the past studies in 
the CCS (from Huyer fit fil** 1989). 
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Figure 2.5 Long-term mean atsKsepheric pressure at sea level 
for (a) January and (b) July (frc»n Huyer, 1983) . 
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Figure 2.6 Long-term (1854-1972) mean wind stress (dynes 
cm’ ) by one-degree squares, from ship reports (Nelson, 
1977). Contours correspond to constant magnitudes of 0.5, 
1.0, and 1.5 dynes cm'^; areas where the wind stress exceeds 
1.0 dynes cm’^ are shaded (from Huyer, 1983). 
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Figure 2.7 An Illustration of the distributions of key 
properties involved in the teleconnection process during a 
typical El Niho-Southem Oscillation event. (Left) Upper 
troposphere anomalies of geopotential height for the 
Northern Hemisphere winter. (Right) Contours of surface 
barometric pressure, with surface wind anomalies (arrows) in 
the tropics (from Enfield, 1989). 
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III. MODEL DESCRIPTION AND EXPERIMENTAL DESIGN 


A. MODEL DESCRIPTION 

1. Basic Model 

The model used In this research is a slightly 
modified version of that used by Batteen r 1. (1989, 1992) 
which was adapted from Haney's (1974, 1985) closed basin 
model. This limited area EBC model is multi-level, uses 
non-adicibatic primitive equations on a jS plane, cuid has both 
barotropic cuid baroclinic coo^onents. The following 
approximations are made: hydrostatic, Boussinesq, euid rigid 
lid. (For a more detailed summary of the model eind the 
governing equations, refer to Batteen st Rl* (1989)). 

2. Friction 

The type of diffusion used Ccua be importcmt for the 
generation of mesoscale eddies. Biharmonic rather them 
Laplacian lateral heat and momentiim diffusion is used in 
this study. The use of biharmonic lateral diffusion allows 
mesoscale eddy generation via barotropic (horizontal-shear) 
or baroclinic (vertical-shear) instability mechanisms, or 
both, as Hollemd (1978) and Holland and Batteen (1986) have 
shown. The biharmonic coefficients are the Scune (i.e., 2.0 
X 10^"^ cm^s"^) as used in Batteen at al* (1989) . Weak (0.5 
cm^s"^) vertical eddy viscosities emd conductivity are used. 
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A simplified quadratic drag law is used to parameterize 
bottom stress (Weatherly, 1972). 

3. Boundary Conditions 

The domain of the model on the eastern boundary is 
closed and is modeled as a straight, vertical wall. A no¬ 
slip condition on the alongshore velocity is used. The 
kinematic boundary condition of no flow through the eastern 
boundary is imposed on the cross-shore velocity component. 

A modified version of the radiation boundary conditions of 
Camerlengo euid O'Brien (1980) is used for the open 
boundaries to the south, west, emd north. 

At the surface, the rigid lid approximation filters 
out external gravity waves and allows a longer time step. 

On the ocean bottom, no vertical heat flux or vertical 
velocity is allowed. 

4. Domain and Resolution 

The domain of the model is an elbow shaped region 
off California's coast from 35** N to 42.5** N as shown in 
Figure 2.3. This covers an area 1408 km alongshore and 1024 
)an in the cross-shore direction. The horizontal resolution 
is 8 km in the cross-shore direction euid 11 km in the 
alongshore direction. Mesoscale features within the CCS 
should be resolvaOsle with this resolution because typical 
eddy wavelengths are around 100 km (Breaker amd Mooers, 

1986). Bottom topography and variations in the coastline 
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have been omitted to focus on the role played by wind 
forcing in the CCS. 

5. Finite Difference Scheme 

A staggered grid model, using the Arakawa auid Lamb 
(1977) B-scheme is used in the horizontal. Variad^les are 
defined as follows: u, v, and w are velocity components 
eastward, northward and vertical, respectively, t denotes 
time, T is temperature, p is density, p is pressure, emd ^ 
is the streamfunction. The center of each grid point 
defines u and v while T, p, w, and p are defined at the 
comers. The B-scheme is used because it gives a superior 
representation of the geostrophic adjustment process 
(Batteen and Heui, 1981) . 

The model uses ten layers in the vertical, with 
constant z-levels at depths of: 13, 46, 98, 182, 316, 529, 
870, 1416, 2283, and 3656 m. This spacing scheme 
concentrates more layers above the thermocline in the 
dynamically active portion of the oceem, consistent with 
Haney (1974). Variables are staggered in the vertical with 
u, V, p, cuid T defined at each level, while w and p are 
defined between the levels which lie between the ocesm 
surface (w^O) euid at the bottom (z—H, where H is the 
depth). A Matsuno time step followed by ten leapfrog time 
steps is used omd continuously repeated during the model 
run. 
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6. Method of Solution 

The method of solution is straightforward with the 
rigid lid and flat bottom assunption because the vertically 
integrated horizontal velocity is subsequently non- 
divergent. The vertical mean flow can be described by a 
strecunfunction rp which can be predicted from the vorticity 
equation, while the vertical shear currents can be predicted 
after the vertical mean flow is subtracted from the original 
equations. The other variables T, p, w, aind p Ceua be 
explicitly obtained from the thermodynamic energy equation, 
equation of state, continuity equation, and hydrostatic 
relation, respectively. 

7. Initial Conditions 

An e3q}onentlal ten^erature profile with a vertical 
length scale of h-450 m is used as the initial meein 
stratification. The bottom tenperature is assumed to be 2^ 

C and the chcuige in temperature between the bottom emd 
surface is 13** C. This profile was derived by Blumberg and 
Mellor (1987) and has been used by Batteen (1989) and 
Batteen £t al- (1989). 

B. WZMD AMD AIMOSPEERZC FORCING 

The winds used to force the model are from the Europecm 
Centre for Medium-Range Weather Forecasts (ECMWF) (Trenberth 
Sii SX-, 1990) cuid have been provided by the National Center 
for Atmospheric Research (NCAR). Winds were extracted from 
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world wide mean monthly wind velocities that were provided 
at 2.5 degree spacing for the years 1980 to 1989 along with 
120 month climatology for the area of study. This data was 
linearly Interpolated In space amd then In time to provide 
dally forcing at the desired grid points. 

The total heat flux across the surface la Initially set 
to zero by choosing an Initial air temperature that forces 
the net flux of longwave radiation, sensible heat, and 
latent heat to zero, as In Batteen al* (1989). The air 
temperature Is then used In the model. Any subsequent 
surface heat flux forcing Is thus a secondary effect of the 
chcinges to SST due to the wind and/or thermal structure of 
the oceaui. 

C. SPECIFIC EZPERIMEIITAL DESIGN 

Climatological winds were used for 365 days (Experiment 
1} and then used to restart the model to run climatological 
winds for another 365 days (Experiment 4). Next, 1981 winds 
were run for 365 days (Experiment 2) and used to Initialize 
the run of 365 days with 1982 winds (Esqperiment 5) . Then 
1983 winds from a cold restart were run for 365 days 
(Esperlment 3) emd then used to initialize the run of 366 
days (leap year) with 1984 winds (Esqperlment 6) . Two 
additional extended runs were conducted in which the model 
rcui for four years. Experiment 7 used the climatological 
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winds four times, while Bscperiment 8 used the monthly winds 
of 1981 through 1984. 

D. BNERGT AKALTSZS TSCHNIQUSS 

The energy technique used is the same as that used auad 
described in Batteen (1992), and is based on that of 

Han (1975) euid Semtner and Mints (1977). This is done to 
gain a better understanding of the flow within the CCS aind 
of the types of energy transfer during iinstable flow. A 
brief summary of the method follows. 

Kinetic energy (K) is calculated for the horizontal 
conponents. After quasi-steady state is reached where the 
total kinetic energy is nearly constant, mean and eddy 
kinetic energies are calculated using the averaged sum of 
squared mecui and eddy horizontal fields, respectively. 

Next, the available potential energy is calculated aad used 
to determine when a quasi-steady state is reached and when 
statistics should be collected. Then, both mean emd eddy 
availcUsle potential energy are conqputed. 

The energy treuasfers follow that of Semtner emd Mlntz 
(1977), and are used to argue for the type of insted)llity 
mechanism (e.g., barotropic, baroclinic, or mixed) which 
leads to the Initial eddy generation in each experiment. 
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ZV. RESULTS OF BZPBRIllBIITS 


Experiments 1 through 3 study the effects of forcing a 
model ocean domain that is initially at rest. E3q>eriment 1 
uses monthly climatological average (1980 - 1989) wind 
fields, while Experiments 2 and 3 use monthly average wind 
fields for 1981 and 1983, respectively. E3q>eriments 4 
through 6 study the effects of forcing a model domain that 
has been spun up by the wind field from the previous year. 
Experiment 4 uses the results of E3q>eriment 1 as the start 
for a second year of climatological average wind fields, 
while Experiment 5 and 6 use the results of Experiment 2 and 
3 as the start for 1982 amd 1984, respectively. Experiments 
7 and 8 use the results of longer (4-year) time runs from 
climatological amd 1981 - 1984 wind fields, respectively, to 
examine the maintenance of features after initial spin up. 

A. EXPERIMENTS STARTED FROM REST 

1. Experlaient 1 - Cllinatologlcal Minds 

Experiment l was forced with monthly averaged 
climatological (1980 - 1989) winds. Figure 4.1 depicts the 
monthly wind patterns over the model domain starting with 
Jamuary 15 (day 15) amd ending on December 15 (day 345). 

The atmospheric pressure pattern for Jamuary has a low to 
the north amd a high to the south which results in a wind 
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divergence near 40® N (Figure 4.1a). This pattern of 
poleward winds north of 40® N and eguatorward winds to the 
south, continues through February and March (Figures 4.1b 
and c) . During April and May (Figures 4.Id auid e), the 
divergence in the wind flow is observed to migrate pole%irard. 
By June (Figure 4.If), an equatorward conponent in the wind 
field is observed along the entire model domain. July, 
August, euid September (Figures 4.1g through i) depict the 
strongest and most equatornirard winds. By October (Figure 
4.1j), the winds start to weaken throughout the domain amd a 
divergence in the wind flow is observed in the north. This 
divergent wind flow pattern continues through November 
(Figure 4.1k). By December (Figure 4.11) the wind 
divergence has returned to around 40® N. 

In response to the prevailing poleward winds in the 
north, a relatively weak poleward coastal surface current of 
~3 cm s'^ was observed by model day 30 (Figure 4.2). It 
retreated further northward \mtil by day 60 (not shown) it 
was not discernible. An eguatorward surface current was 
observed along most of the coast by model day 90 (Figure 
4.3) due to the prevailing eguatorward winds over most of 
the model domain. A poleward xindercurrent, caused by the 
meridional variation of the wind field (Batteen et al.. 
1989), first became evident at model day 60 (Figure 4.4), 
and was fully ested^lished by model day 90 (Figure 4.5). 
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Along with the equatorward surface current, 
upwelling also occurred and was discernible by day 120 
(Figure 4.6). Meanders in the equatorward surface current 
formed by model day 210 amd are evident in both the velocity 
and ten^erature contour plots (Figures 4.7 amd 4.8, 
respectively). By model day 240, meamders began to form 
cold core, cyclonic eddies as seen in the tenqperature 
contour and velocity plots (Figures 4.9 and 4.10, 
respectively) due to vertical and/or horizontal shear 
instad>ilities between the equatorward jet amd the poleward 
undercurrent. The eddies initially appear to be ad>out 30 km 
(Rossby radius of deformation) in diameter auid extend ad)Out 
100 km off the coast, and appear during the period of 
maximum wind stress (July - September). 

The coldest tenperatures due to upwelling were 
reached by model day 240. After model day 270 no upwelling 
was discernible (not shown). At the end of the year, 
meamders amd eddies in a weadc equatorward surface current 
were still present (Figures 4.11 and 4.12). A very weak 
poleward coastal surface current (not shown), amd a poleward 
\indercurrent were also present at the end of the year 
(Figure 4.13). Note that no filaments developed during this 
first climatological year. 

Except for the lack of development of filaments, 
these model results are consistent with availad>le 
observations cited in the literature. For example, the 
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eddies have rotational velocities of around 50-80 cm s*^ 
which Is In agreement with observations of Mooers auid 
Robinson (1984). The poleward undercurrent Is relatively 
weak (~5 cm s*^ or less) and Is In agreement with available 
observations by Huyer fit (1989). 

A time series of total kinetic energy per unit mass 
over the entire domain Is shown In Figure 4.14. In this 
figure there are around three periods when the kinetic 
energy becomes quasi-steady. These time frcunes correspond 
to model days: 207 - 216, 228 - 231, and 297 - 303. Since 
meanders and/or eddies appeared during these time periods, 
more detailed analysis was performed on these Intervals to 
determine the type of Instability mechanism that could 
generate these features. Barotroplc Instaiblllty can result 
from horizontal sh'iar while barocllnlc Instability can 
result from vertical suear In the current. As Figure 4.12 
shows, there Is conslderedsle horizontal cmd vertical shear 
between the equatorward surface current and the poleward 
undercurrent. As a result, both types of Instabilities 
(mixed), can be present simultaneously. Energy transfer 
calculations which consist of barotroplc (mean kinetic 
energy to eddy kinetic energy) emd barocllnlc (eddy 
potential energy to eddy kinetic energy) components were 
performed on the three time frcunes above. 

The results for the lnstcd>illty euialyses for these 
periods are siimmarlzed In Table 4.1 (see i^pendlx A for the 
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energy transfer plots for this experiment), and show the 
following: For model days 207 to 216, 228 to 231, and 297 
to 303 both barotropic and baroclinic (i.e., mixed) 
instabilities were present in the coastal, equatorward 
region of the domain. Additional analysis of model days 297 
to 303 shows barotropic instcd^llities to be weaik (Figure 
A.5) and baroclinic instabilities to be dominant (Figure 
A.6) during this period. 

2. Eaq>erlment 2 - 1981 Winds 

Experiment 2 was forced with the monthly EC34WF winds 
from 1981. Figure 4.15 depicts the monthly wind patterns 
within the model domain starting with Jemuary and running 
through December. Unlike the climatological winds, January 
starts off with poleward flow throughout the model domain 
(Figure 4.15a). February (Figure 4.15b) shows the typical 
wind divergence at 40” N, consistent with the climatological 
data. March through June (Figures 4.15c through f) depict 
the wind divergence moving poleward suid also shows the wind 
speed increasing to a meocimum in June, especially in the 
southern half of the region. Figures 4.15g and h (July and 
August) depict strong winds and equatorward flow throughout 
the model domain. By September (Figure 4.151), the wind 
divergence reappears in the north and migrates equator%^rd 
over the remaining three months of the year (Figures 4.15j 
to 1) . By Deceoiber the wind divergence is at 37.5® N 
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(Figure 4.151), which is around 2.5** farther eguatorward 
than the climatological wind divergence. 

As in Experiment 1, a relatively weak poleward 
coastal surface flow ~2 cm s'^ was observed in the poleward 
region of the model domain (Figure 4.16). In response to 
the wind pattern, an equatorward surface current was 
observed along most of the coast by model day 90 (Figure 
4.17). A weak poleward undercurrent became evident by model 
day 90 (Figure 4.18), and became as large as 5 cm s'^ in 
layer 5 (316 m) by day 210 (Figure 4.19). 

Upwelling of tenperatures cooler them 15** C was 
observed by model day 120 (as observed in the climatological 
run), as seen in Figure 4.20. Mecmders in the surface 
current also formed by the same model day as in the 
climatological run (day 210) as depicted in the velocity amd 
tenperature contour plots (Figures 4.21 and 4.22, 
respectively). By model day 270, the first sign of eddies 
appear (see the velocity amd tenperature contour plots of 
Figures 4.23 amd 4.24), which is ~30 days later tham in the 
climatological run. The amplitude of the meamder of the 
coastal current increases to the south as a result of the 
stronger eguatorward winds in this region. Again, these 
eddies are cold core, cyclonic, amd about 30 km in diameter. 

Upwelling of the coldest surface water oc^'urred by 
model day 240 (not shown), corresponding to the period of 
strongest ecpiatorward winds. Upwelling did not disappear 


30 







until model day 330 (not shown), which is ~60 days later 
than the climatological results. At the end of the year, 
meanders in the equatorward surface current along with 
embedded eddies remain (Figure 4.25). A poleward 
undercurrent is also discernible at the end of the year 
(Figure 4.26) . A very wecdc pole%«ard coastal surface current 
(not shown) has reappeared in the poleward region of the 
model domain. As in Experiment 1, no filaments developed 
during this experiment. 

Figure 4.27 shows the time series of total kinetic 
energy per unit mass over the entire domain for 1981 and 
shows four periods when the kinetic energy becomes quasi- 
steady, model days: 210 - 222, 240 - 252, 297 - 303, euid 
309 - 333. Energy tramsfer plots were obtained for these 
periods amd yielded mixed (i.e., both barotropic amd 
baroclinic) instad>ilities in all cases in the coastal, 
equatorward region of the model domain. These results are 
s\nnmarized in Tadsle 4.1 (see Appendix B for the energy 
transfer plots for this experiment). During the time frame 
297 to 303, barotropic instability (Figure B.5) was much 
greater than baroclinic (Figure B.6) instaibility. 

3. E3q>erlment 3 * 1983 Winds 

Experiment 3 represents the El Nlfio year of 1983 and 
was forced with the monthly winds depicted in Figure 4.28 
representing Jamuary through December. The year starts off 
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with predominantly poleward winds for the first three months 
as seen in Figures 4.28a - c. April shows a change in the 
wind pattern as equatorward winds are present throughout 
(Figure 4.28d). These winds remain equatorward and 
strengthen from May through September (Figures 4.28e • i). 
The winds weaken in October (Figure 4.28j). A divergence in 
the wind field appears near 38** N during November (Figure 
4.28k) and remains through December (Figure 4.281). 

In response to the three consistent months of strong 
poleward winds at the start of 1983, a poleward coastal 
surface current with velocities as high as 5 cm s'^ is 
established by model day 30 (Figure 4.29) and remains 
through model day 60 (not shown). By model day 120 cm 
equatorward surface current develops along most of the coast 
(Figure 4.30). A weak poleward iindercurrent is present at 
model day 90 as seen in Figure 4.31, and strengthens to 
around 6 cm s"^ by model day 195 (not shown) . 

Upwelling and meanders both appear by model day 150 
(Figure 4.32). The meanders are also discernible in the 
velocity plot (Figure 4.33). Note that the meanders appear 
farther poleward than in the previous e3q>eriments. Eddies 
start to form by model day 240 as seen in the velocity 
(Figure 4.34) amd tenperature (Figure 4.35) fields. As in 
the previous eoqperiments, these eddies are cold core and 
cyclonic. 









The coldest upwelllng ten^eratures were reached by 
model day 240 but were warmer (only 14.0** C as conpared to 
13.5* C and 11.5* C for climatology euid 1981, respectively) 
than observed in previous runs. No upwelling was 
discernible after model day 270. By the end of the year, 
meanders in the eguatoirward surface current along with 
embedded eddies (Figure 4.36), a wealc poleward coastal 
current (not shown) , emd a polevrard undercurrent (Figure 
4.37) were discernible. As in the previous e3q>eriments, no 
filaments developed during the eoqperiment. 

Figures 4.38a auid b show the time series of total 
Icinetic energy per unit mass for 1983. Quasi-steady periods 
that pron^ted additional evaluation were, model days: 186 - 
192, 222 - 234, 249 - 258, 279 - 288, 312 - 315, 324 - 327, 
and 342 - 348. The results for the instability analysis for 
this year are summarized in TeUsle 4.1 (see Appendix C for 
the energy trcuisfer plots) . Every period showed mixed (both 
barotropic cuid baroclinic) inst8d}ilities to be present. The 
instability analysis for model day 279 to 288 was the only 
case where one type of instability was noticeably greater 
than the other, and in this case barotropic (Figure C.7) was 
greater them baroclinic (Figure C.8). 

4. Summary of First Tsar Run Results 

Ted}le 4.2 summarizes the climatological, 1981, and 
1983 experiment runs for the first year. In all first year 
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runs, both a surface equator%«ard current and a poleward 
undercurrent developed. There was also a poleward coastal 
surface current present during the winter season. Upwelling 
was a seasonal feature which, as expected, had the coolest 
temperatures during a non-El Nifio year (1981) and warmest 
tenperatures during an El Nlho year (1983) . Both mecuiders 
and eddies were generated. The eddies were cold core, 
cyclonic, and had a size aroimd the Rossby radius of 
deformation (30 km). Along with an equatorward surface 
current, both mecmders amd eddies remained at the end of 
each year. Instability auialyses revealed that both 
barotropic and baroclinic processes were inportouit eddy 
generatiox. mechcmisms. In all first year inins, no filciments 
were observed. 

B. EXPERIMENTS SPUN UP FROM TEE PREVIOUS TEAR 
1. Experiment 4 - Climatological Ninds 

Experiment 4 was forced with the same monthly 
averaged climatological winds as in Experiment 1 (Figure 
4.1). In this experiment, the model was continued with the 
spun-up output from Experiment 1, rather than starting from 
rest. For convenience this experiment is referred to as 
"Clirno II" because climatological winds are used a second 
time. 

In response to the wind forcing, an eguatoirward 
surface current was observed by model day 396 (year day (yd) 
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30) (Figure 4.39) . Entbedded In this current are meanders 
and cyclonic eddies (Figures 4.39 and 4.40). A weak 
poleward coastal surface current was also present initially 
(not shown) . A poleward undercurrent was estcdslished by 
model day 426 (yd 60) (Figure 4.41). 

Upwelling was observed by model day 546 (yd 180) 
(Figure 4.42), and peaked by model day 606 (yd 240) (Figure 
4.43), with a minimum ten^erature near 12C. Cessation of 
upwelling occurred by model day 726 (yd 360) (not shown). 

In previous model runs anticyclonic mecinders were 
observed, but they did not close off to form emticyclonic 
eddies. During this experiment, in addition to the 
predominantly cyclonic eddies, a warm core euiticyclonic eddy 
(defined by closed contour in streamline plot - not shown) 
formed by model day 456 (yd 90) (Figure 4.44). This 
anticyclonic eddy persisted for the remainder of the year. 

Figure 4.45 is a cross section in the middle of the 
domain and depicts the broad meanders in the CC along with a 
strong core of the CUC near 250 m in depth. Between model 
days 696 (yd 330) euid 726 (yd 360), an eddy formed as the 
result of a pinched-off mecinder (Figures 4.46 eUid 4.47). A 
fllcunent was observed during this experiment by day 606 (yd 
240) oind extended about 300 km offshore (15* C isotherm) 
near 40* N (Figure 4.43) . 

A time series of total kinetic energy per unit mass 
over the entire domain is shown in Figure 4.48. Four 
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perlcxis were investigated further because the kinetic enex^ 
became quasi-steady. These time frames correspond to model 
days: 378 - 408 (yd 12 - 42), 453 - 465 (yd 87 - 99), 504 - 

516 (yd 138 - 150), and 624 - 636 (yd 258 - 270). All four 
of these periods consisted of both barotropic and baroclinic 
instabilities. These results are summarized in TcJsle 4.3 
(see Appendix D for the energy trauisfer plots for this 
e3q>eriment) . 

2. Experiment 5 - 1982 Winds 

Experiment 5 was forced with monthly winds from 
1982. This experiment continues the spun up case using the 
output of Experiment 2 (1981) as a starting point, rather 
theui starting from rest. Figure 4.49 shows the monthly wind 
patterns within the model domain starting with January and 
running through December. January starts with the typical 
divergent flow, this year slightly to the north at 42® N 
(Figure 4.49a). February moves the divergence eguatorward 
to 40® N and has weedcer winds to the south (Figure 4.49b) . 
The flow in March is onshore throughout the model domain 
with a slight poleward conponent in the north (Figure 
4.49c). April has the divergence at 42.5® N (Figure 4.49d). 
May starts the all eguatorward flow and is the strongest 
(Figure 4.49e), this regime continues in J\me (Figure 4.49f) 
through October (Figure 4.49i). A divergent flow returns in 
October at 42.5® N (Figure 4.49j), euid migrates eguatorward 
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to 38.5** N in November (Figure 4.49k) and remains in 
December (Figure 4.491). 

As in the previous experiment, the year starts with 
a fully developed equatorward surface current which contains 
embedded mecuiders and eddies (Figure 4.50). Also present is 
an anticyclonlc eddy, euid a poleward surface coastal current 
in the northern part of the model domain (Figures 4.50 and 
4.51). A relatively strong (~9.5 cm s*^) poleward 
undercurrent is present by model day 426 (yd 60). A cross 
section in the middle of the domain (Figure 4.52) shows this 
feature along with a strong subsurface flow associated with 
an eddy (Figure 4.53). 

Upwelling begain by model day 516 (yd 150) (Figure 
4.54), euid pecUced by model day 606 (yd 240) (Figure 4.55). 
The coolest water observed was around 12.0® C. The 
upwelling ceased by model day 696 (yd 330) (Figure 4.56). 

Four cyclonic eddies eUid one cmticyclonic eddy were 
presei: (e.g.. Figure 4.53) throughout most of this 

experiment. The size ranged between 30 cUid 150 km in 
diameter and the westward propagation of these larger 
features was aroiind 500 km. As in Experiment 4, a filament 
was observed during this experiment. By day 600 (yd 234) 
(not shown) this fllcunent extended about 350 km offshore 
near 38.5** N. A cold core eddy pinched off from this 
filament by day 606 (yd 240) (Figure 4.55). 
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A time series of total kinetic energy per unit mass 
over the entire domain was produced (Figure 4.57), and 
showed two quasi-steady periods for additional study. These 
were model days: 507 - 516 (yd 141 - 150), and 525 - 528 
(yd 159 - 162). As in the previous experiments, the energy 
transfer plots indicate the presence of both barotroplc and 
baroclinic instcUOilities. Table 4.3 summarizes these 
results (see J^pendix E for the energy transfer plots for 
this experiment). 

3. E3q>eriment 6 - 1984 Winds 

Experiment 6 was forced with monthly winds from 
1984. This experiment continues the spun-up case using the 
output of Experiment 3 (1983) as a starting point, rather 
than starting from rest. Figure 4.58 shows the monthly wind 
patterns within the model domain starting with January amd 
running through December. January has a divergent flow at 
43** N (Figure 4.58a), while February shows that the 
divergence has moved equatorward to 39* N (Figure 4.58b). 
From March through May (Figures 4.58c through e), the wind 
divergence moves from 41* N to 43* N and the winds become 
stronger. June is the start of equatorward wind flow 
(Figure 4.58f), which continues through July and August 
(Figures 4.58g aind h) . In September, the wind divergence 
reappears in the north at 47* N (Figure 4.581). During 
October auid November, the divergence migrates equatorward to 
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37** N (Figures 4.58i and k) . The year ends in December with 
weak winds along the coast with a divergence at 40^ N, amd 
with predomincuitly eguatorward winds further off the coast 
(Figure 4.581). 

At the start of the experiment, there is a well- 
developed eguatorward surface flow with embedded meanders 
and eddies (Figure 4.59). There is also a weaUc poletfard 
coastal surface current in the northern comer of the model 
domain which remains unti] model day 426 (yd 60) (Figure 
4.60). A relatively strong (**10 cm s‘^) poleward 
undercurrent is evident, as seen in a cross-section in the 
middle of the model domain, at model day 546 (yd 180) 

(Figure 4.61). 

Upwelling begins by model day 456 (yd 90) (Figure 

4.62) . Corresponding to the period of maximum eguatorward 
winds, the peak occurs at model day 576 (yd 210) (Figure 

4.63) in which cooler ten^eratures of around 10** C are 
present along the entire coast. This year had upwelling of 
cooler water earlier thoui previous years emd also had the 
coldest water which was a result of the strong eguatorward 
winds. Cessation of the cooler water occurs by model day 
726 (yd 360) (not shown). 

The eddies are around 90 km in dieuneter and stay 
within 300 km of the coast. A filament occurs by model day 
666 (yd 300) (Figure 4.64). This filament is observed near 
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38** N as a plume of 15** C water aid extends about 200 km off 
the coast. 

The time series plot of total kinetic energy per 
unit mass over the entire domain Is shown In Figure 4.65. 

The kinetic energy Is still fluctuating because the model Is 
still spinning up during the first part of the year. Three 
periods were selected for additional study, model days: 576 
- 591 (yd 210 - 225), 606 - 624 (yd 240 - 258), and 645 - 
657 (yd 279 - 291). All of these Included both barotroplc 
cuid barocllnlc Instabilities, and are summarized In Tcible 
4.3 (see Appendix F for the energy treuisfer plots for this 
experiment) . 

4. Sunmazy of Second Tear Rim Results 

Table 4.4 summarizes the follow-on (second year) 
model runs. In general, a weaJc coastal poleward surface 
current was evident during the winter season In the northern 
region of the domain for all years (Cllmo II, 1982 euid 
1984) . A meeuiderlng eguatorvmrd surface current, already 
present, at the start of these esqperlments was maintained. 

A poleward imdercurrent was generated In all experiments. 

The eddies, already present at the start of the experiments, 
were both cold core, cyclonic, euid warm core, cuiticyclonlc 
eddies In Cllmo II euid 1982, while in 1984 there were only 
cold core, cyclonic eddies. The size and offshore extent of 
the eddies were consistent with the first e3q)erlmental runs. 
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Seasonal upwelllng occurred with minimum temperatures of 
~12® C for Cllmo II euid 1982, and of -10® C for 1984. The 
year 1984 had the largest extent of cool upwelled water as 
expected for a year with the strongest eguatorward winds. 
Filaments developed in all three esqperiments, but was most 
pronounced (as expected) in 1984. The filaments had cui 
offshore extent between 200 eind 350 km. Insted>ility 
analyses revealed that both barotropic amd baroclinlc 
processes were involved in eddy generation. 

5. Bacperlment 7 - Multiple Climatological Runs 

Experiment 7 was aui extended climatological run in 
which the climatological winds used in E3q>eriment 1 (Figure 
4.1) were used in repetition four times. The third auid 
fourth years of the run follow Bacperiment 4 (Climo II) amd 
are dubbed "Climo III" and "Climo IV", respectively. (The 
values in the model on the last day of a run were retained 
for use at the start of the next run, rather tham starting 
from scratch or rest as is done in a cold restart.) 

In response to the wind forcing, am eguatorward 
surface current was observed by model day 760 amd 1125 (yd 
30) for Climo III amd Climo IV, respectively (not shown). 
Embedded in this current were meanders amd cyclonic eddies 
(not shown) . A weaik poleward coastal surface current was 
also present initially (not shown). A poleward imdercurrent 
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was evident early on by model day 820 and 1185 (yd 60) for 
Climo III emd Climo TV, respectively (not shown). 

Upwelling was observed by model day 850 eind 1215 (yd 
120) for Climo III and Climo IV, respectively (not shown). 
The upwelling peaiked by model day 970 (yd 240) for Climo III 
(not shown) and day 1305 (yd 210) for Climo IV (not shown), 
with a minimum tenperature near 11** C. Cessation of 
upwelling occurred by model day 1060 (yd 330) for Climo III 
(not shown), and was still present in small patches by day 
1455 (yd 360) for Climo TV (not shown). 

Figure 4.66 and 4.67 show the total kinetic energy 
per unit mass time series over the entire domain for Climo 
III and Climo IV, respectively. Climo III has a peak 
kinetic energy of 166 cm^s’^ around model day 915 (yd 185) 
while Climo IV has a peak of 145 cm^s*^ at model day 1293 
(yd 198). Figure 4.68 depicts the total kinetic energy per 
unit mass time series for all four climatological runs 
together. It clearly shows that the first emd second runs 
were involved in the generation or spin-up of features in 
the model, while the third and fourth runs oscillated eibout 
125 cm^s'^ and indicate the maintenance of these features 
(i.e., a quasi-steady state). 

6. Experiment 8 - Multiple Tears (1981 - 1984) 

Experiment 8 was a run of two additional years of 
the model using the monthly winds of 1983 (Figure 4.28) euid 
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1984 (Figure 4.58). This experiment differs from Experiment 
3 (1983) which was started from scratch (cold restart). The 
runs in this experiment followed in succession after 
Experiment 5 (1982), cuid used the values in the model on the 
last day of the year as input to the restart. 

In response to the wind forcing, etn eq^aatorward 
surface current was observed by model day 760 auid 1125 (yd 
30) for 1983 and 1984, respectively (not shown). Embedded 
in this current were meanders and cyclonic eddies (not 
shown). A weak poleward coastal surface current was also 
present initially (not shown). A poleward undercurrent was 
evident by model day 820 cuid 1185 (yd 60) for 1983 amd 1984, 
respectively (not shown). 

Upwelling was observed by model day 880 (yd 150) for 
1983 (not shown), and day 1155 (yd 60) for 1984 (not shown). 
The upwelling peaiked by model day 940 auid 1305 (yd 210) for 
1983 and 1984, respectively (not shown). The minimum 
ten^erature was near ll^* C for both 1983 amd 1984, but 1983 
had only patches of this water while 1984 had this cool 
water along the entire coast. Upwelling was still present 
in small patches by day 1090 amd 1455 (yd 360) for 1983 amd 
1984, respectively (not shown). 

Figure 4.69 amd 4.70 are the total kinetic energy 
per unit mass time series over the entire domain for 1983 
and 1984, respectively. The year 1983 has a peadc kinetic 
energy of 128 cm^s'^ around model day 963 (yd 233) amd 1984 
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has a peak of 215 cm^s'^ at model day 1368 (yd 273). Figure 
4.71 depicts the runs for all four years (1981 through 1984) 
together on one total kinetic energy per unit mass time 
series. It appears that the first year (1981) as seen in 
Experiment 2 and the second year (1982) as seen in 
E3q)eriment 5, were involved in generation or spinning up of 
the model. The third year (1983), which was cm El Nifio 
year, had an energy level on par with 1982, so that no 
increase in kinetic energy occurred. In the fourth year 
(1984), the energy level Increase significemtly to the 
highest value (215 cm^s*^) . Conparing this figure to the 
climatological time series (Figure 4.68) for the last two 
years shows that the oieiximum (1984) emd minimum (1983) 
kinetic energies lie on either side (as expected) of the 
climatological year. During a non-El Nifio year, as is 1984, 
features such as equatorward currents auid upwelling would 
strengthen while during an El Nifio year, as in 1983, these 
features would be weaker. 

7. Summary of Multiple Tear Run Results 

Table 4.5 summarizes the last two years of the 
extended (four year) runs. Listed here are Climo IV, 1983, 
cmd 1984. Similar results are obtained: A poleward coastal 
surface current, a surface equatorward current, a poleward 
undercurrent, along with embedded mecuiders and eddies. Of 
interest, upwelling occurs later in 1983 (at yd 150) them 
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for the other years (yd 60 to yd 90). Although the minimum 
temperature is the scune (11.0® C) for all years, 1904 
experienced a larger area of cooler water (all along the 
coast) than that of the El Nifio (1983) year. In all 
experiments, there were both cold amd warm core eddies. 
Finally, filaments were observed during Climo IV euad 1984, 
but not during 1983. 


45 





TABLB 4.1 SUMMARY OF INSTABILITIES FOR CLIMATOLOGICAL YEAR 
(TOP), 1981 (CENTER), AND 1983 (BOTTOM). 


MODEL DAYS 
for 

Climatology 

FEATURE TYPE 
(Meander, 

Eddy, Both) 

INSTABILITY 

(Barotropic, 

Baroclinic, 

Mixed) 

LOCATION 

(Coastal, 

Offshore, 

Both) 

207-216 

Meander 

Mixed 

Coastal 

228-231 

Meander 

Mixed 

Both 

297-303 

Both 

Mixed 

Both 


MODEL DAYS 
for 

1981 

FEATURE TYPE 
(Meander, 

Eddy, Both) 

INSTABILITY 

(Barotropic, 

Baroclinic, 

Mixed) 

——1 

LOCATION 

(Coastal, 

Offshore, 

Both) 

210-222 

Meander 

Mixed 

Coastal 

240-252 

Meander 

Mixed 

Both 

297-303 

Both 

Mixed 

Both 

309-333 

Both 

Mixed 

Both 


MODEL DAYS 
for 

1983 

FEATURE TYPE 
(Meander, 

Eddy, Both) 

INSTABILITY 

(Barotropic, 

Baroclinic, 

Mixed) 

LOCATION 

(Coastal, 

Offshore, 

Both) 

186-192 

Meander 

Mixed 

Coastal 1 

222-234 

Meander 

Mixed 

Both 1 

249-258 

Meouider 

Mixed 

Coastal 

279-288 

Meamder 

Mixed 

Coastal 

312-315 

Both 

Mixed 

Both 

324-327 

Both 

Mixed 

Both 

342-348 

Both 

Mixed 

Both 
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TABLE 4.2 SUMMARY OF EXPERIMENTS - FIRST YEAR. 


FEATURE 

CURRENTS! _ 

Sfc Poleward 

Sfc Equatorward 

Undercurrent 
(poleward only) 

UPNELLZNGt 
(<15 C) _ 

Start _ 

Mciximum 

(min temp) _ 

Cessation _ 

MEANDERS! _ 

Start _ 

EDDIES! _ 

Start _ 

Core _ 

Rotation _ 

Size (diameter) 
FILAMENTS! 



330-360 days 330-360 days 330-360 days 

~90 days _ ~90 days _ ~120 days 

-60 days -120 days -60 days 



-120 days _ -120 days _ -150 days 

-240 days -240 days -240 days 


(13.0® C) 


(11.0* C) 


(14.0® C) 


-270 days 


-330 days 


-270 days 


-210 days 


-210 days 


-150 days 


-240 days_| -270 days 


cold 

cold 

cyclonic 

cyclonic 

-30-60 km 

-30-60 km 



-240 days 
cold 



Start 

Offshore 

extent 


n/a 


n/a 


n/a 






































































TABLE 4.3 SUMMARY OP INSTABILITIES FOR CLIMATOLOGICAL II YEAR 
(TOP), 1982 (CENTER), AND 1984 (BOTTOM). 


MODEL DAYS 
for 

Climo II 

FEATURE TYPE 
(Meander, 

Eddy, Both) 

INSTABILITY 

(Barotropic, 

Barocllnic, 

Mixed) 

LOCATION 

(Coastal, 

Offshore, 

Both) 

378-408 

(12-42) 

Both 

Mixed 

Both 

453-465 

(87-99) 

Meander 

Mixed 

Coastal 

504-516 

(138-150) 

Both 

wixed 

Both 

624-636 

(258-270) 

Both 

Mixed 

Both 


MODEL DAYS 

FEATURE TYPE 

INSTABILITY 

LOCATION 

for 

(Meander, 

(Barotropic, 

(Coastal, 

1982 

Eddy, Both) 

Baroclinic, 

Offshore, 


Mixed) 

Both) 

507-516 

(141-150) 

Both 

Mixed 

Both 

525-528 

(159-162) 

Both 

Mixed 

Both 1 


MODEL DAYS 
for 

1984 

FEATURE TYPE 
(Meamder, 

Eddy, Both) 

INSTABILITY 

(Barotropic, 

Baroclinic, 

Mixed) 

LOCATION 

(Coastal, 

Offshore, 

Both) 

576-591 

(210-225) 

Meauider 

Mixed 

Coastal 

606-624 

(240-258) 

Meander 

Mixed 

Coastal 

645-657 

(279-291) 

Meander 

Mixed 

Coastal 
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TABLE 4.4 SUMMARY OF EXPERIMENTS - SECOND YEAR. 


I 


FEATURE 

CLIMO II 

1982 

1984 

CURRENTS: 


. 


Sfc Poleward 

-366-426 
(1-60), & 
696-726 
(330-360) 
days 

-366-426 
(1-60), & 
696-726 
(330-360) 
days 

-366-426 
(1-60), & 
696-726 
(330-360) 
days 

Sfc Equatorward 

-396 days 
(30) 

-396 days 
(30) 

-396 days 
(30) 

Undercurrent 
(poleward only) 

-426 days 
(60) 

-426 days 
(60) 

-546 days 
(180) 

UPNELLINGt 
(<15 ® C) 




Start 

-546 days 
(180) 

-516 days 
(150) 

-456 days 
(90) 

Maximxim 
(min temp) 

-606 days 
(240) 

(12.0« C) 

-606 days 
(240) 

(12.0® C) 

-576 days 
(210) 

(10.0® C) 

Cessation 

-726 days 
(360) 

-696 days 
(330) 

-726 days 
(360) 

MEANDERS: 




Start 

-396 days 
(30) 

-396 days 
(30) 

-396 days 
(30) 

EDDIES: 




Start 

-396 days 
(30) 

-396 days 
(30) 

-396 days 
(30) 

Core 

cold i warm 

cold & warm 

cold 

Rotation 

cyclonic & 
anticyclonic 

cyclonic & 
euiticyclonic 

cyclonic 

Size (diameter) 

-30-90 )an 

-30-150 km 

-30-90 km 

FILAMENTS: 




Start 

-606 days 
(240) 

-600 days 
(234) 

-666 days 
(300) 

Offshore 

extent 

-300 km 

-350 Jon 

-200 3an 
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TABLE 4.5 SUMMARY OF EXPERIMENTS - MULTIPLE YEARS. 


FEATURE 

11111^22201111 

1983 

1984 

CURRENTS: 




Sfc Poleward 

-1096-1215 
(1-120), & 
1395-1455 
(300-360) 
days 

-731-850 
(1-120), & 
1060-1090 
(330-360) 
days 

-1096-1185 
(1-90), & 
1425-1455 
(330-360) 
days 

Sfc Eguatorward 

-1125 days 
(30) 

-760 days 
(30) 

-1125 days 
(30) 

Undercurrent 
(poleward only) 

-1155 days 
(60) 

-790 days 
(60) 

-1155 days 
(60) 

UPNELLINGt 
(<15 « C) 




Start 

-1185 days 
(90) 

-880 days 
(150) 

-1155 days 
(60) 

Meiximum 
(min ten^) 

-1305 days 
(210) 

(11.0« C) 

-940 days 
(210) 

(11.0» C) 

-1305 days 
(210) 

(11.0“ C) 

Cessation 

-1455 days 
(360) 

-1090 days 
(360) 

-1455 days 
(360) 

MEANDERS: 

' 



Start 

-1125 days 
(30) 

-760 days 
(30) 

-1125 days 
(30) 

EDDIES: 




Start 

-1125 days 
(30) 

-760 days 
(30) 

-1125 days 
(30) 

Core 

cold & warm 

cold & warm 

cold & warm 

Rotation 

cyclonic & 
anticyclonic 

cyclonic & 
euiticyclonic 

cyclonic & 
anticyclonic 

Size (diameter) 

-30-80 Km 

-30-80 km 

-30-80 km 

FILAMENTS: 




Start 

-1305 days 
(210) 

none 

-1305 days 
(210) 

Offshore 

extent 

-350 km 

n/a 

-350 km 
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Figure 4.1 Mid-monthly averaged climatological (1980 - 
1989) ECMWF winds in m s'^: (a) January, (b) February, (c) 

March, (d) April, (e) May, (f) June, (g) July, (h) August, 
(i) September, (j) October, (k) November, and (1) December. 
The latitudinal grid point 1 (129) corresponds to 35® N 
(47.5® N); so that the latitudinal grid point of 52 
corresponds to 40® N. Maximum wind vector is 20 m s*^. 
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Figure 4.1 (continued). 
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Figure 4.1 (continued). 
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Distance off shore (km) 
-853.3 -^682.7 -512 8 -3-11. 


1488.0 


1173.3 


Velocity at model day 90.0 


nitximiH vrCinn 


Figure 4.3 Experiment 1: Surface velocity vectors at day 
90. To avoid clutter, velocity vectors are plotted at every 
third grid point in both the cross-shore and alongshore 
direction, and velocities less than 5 cm s'^ are not 
plotted. 





Depth (m) 


Distance off shore (km) 



Figure 4,4 Experiment 1; Cross-section at y-704 km 
(41.25® N) of the meridional component of velocity at day 
60. Dashed lines Indicate equatorward flow, while solid 
lines Indicate poleward flow. 







Depth (m) 


Distance off shore (km) 

- 053.3 - 602.7 - 512.0 - 341.3 - 170.7 0.0 



North-South velocity at model day 90.0 


Figure 4.5 Experiment 1: Cross-section at 7-704 km 
(41.25* N) of the meridional consonant of velocity at day 
90. Dashed lines indicate equatorward flow, while solid 
lines indicate poleward flow. 
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Temperature contour at model day 120.0 


Figure 4.6 Experiment 1: Surface temperature at day 120 
The contour Interval Is C. The temperature decreases 
towards the coast. 
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Velocity at model day 210.0 


N»*iMii« vrTinn 


Figure 4.7 Experiment 1: Surface velocity vectors at day 
210. To avoid clutter, velocity vectors are plotted at 
every third grid point in both the cross-shore and 
alongshore direction, and velocities less than 5 cm s'^ are 
not plotted. Note start of meander at y~500 km. 







Distacce alongshore (km) 


Distance off shore (km) 

-n53 3 -603.7 -513.0 -3-11.3 -1 70.7 0.0 



Temperature contour at model day 210.0 


Figure 4.8 Experiment 1: Surface temperature at day 210. 
The contour interval is 1* C. The temperature decreases 
towards the coast. Note start of meander at y~500 km. 
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Temperature contour at model day 2-10.0 


Figure 4.9 Experiment 1: Surface temperature at day 240. 
The contour interval is 1® C. The temperature decreases 
towards the coast. Note the formation of cold core eddy at 
y~230 km. 
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Temperature contour at model day 360.0 


Figure 4.11 Experiment 1: Surface temperature at day 360. 
The contour interval is C- Note meanders and eddies are 
still discernible in temperature contours. 








Distance alongshore (km) 



Figure 4.12 Experiment 1: Surface velocity vectors at day 
360 as in Figure 4.3. Note equatorward surface current over 
a poleward undercurrent. The number associated with the low 
at y~800 km corresponds to ten^erature, not velocity. 
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Norlh-Soulh velocity al model day 360.0 


Figure 4.13 Experiment 1: Cross-eection of the meridional 
component of velocity at day 360. Note meanders and eddies 
are still discernible in a weak equatorward surface current 
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Figure 4.14 Experiment 1: Total kinetic energy per unit 
mass time series for climatological (1980-1989) year over 
the entire domain (Units of the kinetic energy are cm^s'^) 
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Figure 4.15 Mid-monthly ECMWP winds for 1981 in m s'^; 

(a) January, (b) February, (c) March, (d) April, (e) May, 

(f June, (g) July, (h) August, (i) September, (j) October, 
(ic) November, and (1) December. The latitudinal grid point 
1 (129) corresponds to 35° N (47.5° N); so that the 
latitudinal grid point of 52 corresponds to 40° N. Maximum 
wind vector is 20 m s‘^. 
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Figure 4.15 


(continued). 
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Figure 4.15 (continued). 
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Figure 4,17 Experiment 2: Surface velocity vectors at day 
90. To avoid clutter, velocity vectors are plotted at every 
third grid point in both the cross-shore and alongshore 
direction, and velocities less than 5 cm s*^ are not 
plotted. 
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Depth (m) 


Distance off shore (knn) 


-102-1.0 -053.3 -602.7 -512.0 -341.3 -170.7 0.0 



North-South velocity at model day 90.D 


Figure 4.18 Experiment 2; Cress-section at 7-704 km 
(41.25® N) of the meridional consonant of velocity at day 
90. Dashed lines indicate equatorward flow, while solid 
lines indicate poleward flow. Note weak poleward 
undercurrent. 
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Distance alongshore (km) 



Figure 4.19 Experiment 2; Layer 5 (316 m) velocity vectors 
at day 210. The maximum poleward velocity is 5 cm s'^ near 
y~469 km. 


73 




Temperalure conlour al model day 120.0 


Figure 4.20 Experiment 2: Surface temperature at day 120 
The contour Interval Is !• c. The temperature decreases 
towards the coast. 












Distance alongshore (km) 


ie?4 e 


Men e 


-B53 3 

_l_ 


Distance off shore (km) 
-682 7 - 512 # -341 3 


1 1 73 3 


/ 

/ / 
/ / 
/ / / 
'II 
'II 

/ / 

'II 

//'' 

' ! / 

' //' 

''ll 

' y/y-' 

' ' ' ! 


! / 


'/ 


/ 


23-1 7 


UJ/^' fll 

/////✓ /I 

-'////^ I 

///// / 

///// 

^///// 

nut 


Velocity at model day 210 0 


I 


tedtlMW vrfifi 


Flgur* 4.21 Experiment 2: Surface velocity vectora at day 
210. Note the meander near y^dSS km. 
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Temperalure contour at model day 210.0 


Figur* 4.22 Experiment 2: Surface ten^erature at day 210 
The contour interval Is 1* C. The ten^erature decreases 
towards the coast. Note the meander near y~469 km. 
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Figure 4.23 Experiment 2: Surface velocity vectors at day 
270. The number associated with the high at y~150 km 
corresponds to tenqperature, not velocity. Note the eddy 
near y~600 km. 
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Temperature contour at model dey 270.0 


Figure 4.24 Experiment 2: Surface temperature at day 270. 
The contour interval Is 1® C. The temperature decreases 
towards the coast. Note the cold core eddy near y~6000 km. 
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Tempernlure contour at model 6ity 360.0 


Flgura 4.25 Experiment 2: Surface temperature at day 360. 
The contour interval la !• C. The temperature decreases 
towards the coast. Note existence of meanders and embedded 
eddies. 
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Flgur* 4.2S Ebcperiment 2: Cross-section at y~704 km 
(41.25* N) of the meridional component of velocity at day 
360. Dashed lines indicate equatorward flow, while solid 
lines indicate polemrd flow. Note existence of poleward 
undercurrent. 


80 









{ 3*6/2*,U.3 . 3> 


Kirietic Energy vs. Time 



Flgur* 4.27 Experiment 2: Total kinetic energy per unit 
mass time series for 1981 over the entire domain (Units of 
the kinetic energy are cm^s*^) . Quasi-steady periods: days 
210-222, 240-252, 297-303, and 309-333. 
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Figure 4.28 (continued). 
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Figure 4.28 (continued). 
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Figure 4.30 Experiment 3: Surface velocity vectors at day 
120. To avoid clutter, velocity vectors are plotted at 
every third grid point In both the cross-shore and 
alongshore direction, and velocities less than 5 cm s'^ are 
not plotted. Note equatorward surface current dominates. 

The number associated with the high at y~250 km corresponds 
to temperature, not velocity. 
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Distance off shore (km) 
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North-South velocity at model day 90.0 


Figure 4.31 Experiment 3: Cross-section at y-704 km 
(41.25® N) of the meridional component of velocity at day 
90. Dashed lines indicate eguatorward flow, while solid 
lines indicate poleward flow. 
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Temperature contour at model day 150.0 


Figure 4.32 Experiment 3: Surface tenperature at day 150 
The contour interval is 1® c. The ten^erature decreases 
towards the coast. Note the existence of upwelling and 
meanders. 
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Figure 4.33 ETcperlment 3: Surface velocity vectors at day 
150. To avoid clutter, velocity vectors are plotted at 
every third grid point in both the cross-shore and 
alongshore direction, and velocities less than 5 cm s'^ are 
not plotted. The number associated with the high at y~850 
km corresponds to temperature, not velocity. Note the 
meanders near y~938 km. 
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Figure 4.34 Experiment 3: Surface velocity vectors at day 
240. To avoid clutter, velocity vectors are plotted at 
every third grid point in both the cross-shore and 
alongshore direction, and velocities less than 5 cm s'^ are 
not plotted. The number associated with the high at y~250 
km corresponds to temperature, not velocity. Note existence 
of eddies. 










Temperature contour at model day 240.0 


Flgura 4.35 Experiment 3: Surface temperature at day 240 
The contour interval is 1* C. The temperature decreases 
towards the coast. Note existence of eddies. 
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Temperature contour el model day 360.0 


Flgur« 4.3€ Bsqperlment 3: Surface temperature at day 360. 
The contour interval Is 1** C. The tenperature decreases 
towards the coast. Note eddies remain at end of year. 
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Flgur* 4.37 Experiment 3; Cross-section at y-704 km 
(41.25* N) of the meridional con^nent of velocity at day 
360. Dashed lines indicate equatorward flow, while solid 
lines indicate poleward flow. Note poleward undercurrent 
remains at end of year. 
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Kinetic Energy vs. Time 



TI m» I days ) 


Figure 4.38a Experiment 3: Total kinetic energy per unit 
mass time series for 1903 (days 1-306) over the entire 
domain (Units of the kinetic energy are cm^s'^) . Quasi- 
steady periods: days 186-192, 222-234, 249-258, and 279-288. 
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Figure 4.39 Experiment 4: Surface velocity vectors at day 
396 (yd 30). To avoid clutter, velocity vectors are plotted 
at every third grid point in both the cross-shore and 
alongshore direction, and velocities less than 5 cm s'^ are 
not plotted. The number associated with the low at y~0O km 
corresponds to teit^ierature, not velocity. Note equatorward 
surface current with embedded meanders and cyclonic eddies. 
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Temperature contour at mode) day 390.0 


Figure 4.40 Experiment 4: Surface teziperature at day 396 
(yd 30). The contour interval is 1® C. The temperature 
decreases towards the coast. Note cold core eddies. 
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Figure 4.41 E3q>erlinent 4: Cross-section at y-704 km 
(41.25* N) of the meridional ccxi^onent of velocity at day 
426 (yd 60). Dashed lines indicate eguatotvrard flow, while 
solid lines indicate poleward flow. Note existence of 
poleward undercurrent. 
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Temperalure contour at model day 546.0 


Figure 4.42 Experiment 4: Surface ten^erature at day 546 
(yd 180). The contour interval is 1« C. The tenperature 
decreases towards the coast. 
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Temperature contour at model day 606.0 


Figure 4.43 Experiment 4: Surface temperature at day 606 
(yd 240). The contour interval is 1® C. The temperature 
decreases towards the coast. 















Distance alongshore (km) 


Distance off shore (km) 



vrfion 


Figure 4.44 Experiment 4: Surface velocity vectors at day 
456 {yd 90). To avoid clutter, velocity vectors are plotted 
at every third grid point in both the cross-shore and 
alongshore direction, and velocities less than 5 cm s'^ are 
not plotted. Note anticyclonic eddy at y~650 km. 
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North-South velocity at model day 666.0 


Figure 4.45 Experiment 4 : Cross-section at y-704 km 
(41.25® N) of the meridional component of velocity at 
666 (yd 300). Dashed lines indicate eguatorward flow 
solid lines indicate poleward flow. 
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Temperature contour at mode! day 606.0 


Figure 4.46 Experiment 4: Surface ten^erature at day 696 
(yd 330), The contour Interval is 1® C. The temperature 
decreases towards the coast. Note formation of eddy from 
pinched off meander at y~469 km. 
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Figure 4.48 Experiment 4: Total kinetic energy per unit 
mass time series for Climo II over the entire domain (Units 
of the kinetic energy are cm^s'^) . Quasi-steady periods: 
days 378-400 (yd 12-42), 453-465 (yd 87-89), 504-516 (yd 
138-150), and 624-636 (yd 258-270). 
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Figure 4.49 Mid-monthly EQWF winds for 1982 in m s'^: 

(a) January, (b) February, (c) March, (d) April, (e) May, 

(f) June, (g) July, (h) August, (i) SeptenOser, (j) October, 
(k) November, and (1) December. The latitudinal grid point 
1 (129) corresponds to 35** N (47.5** N); so that the 
latitudinal grid point of 52 corresponds to 40** N. Maximum 
wind vector is 20 m s'^. 
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Figure 4.49 (continued) 
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Figure 4.50 Experiment 5: Surface velocity vectors at day 
396 (yd 30). To avoid clutter, velocity vectors are plotted 
at every third grid point in both the cross-shore and 
alongshore direction, and velocities less than 5 cm s‘^ are 
not plotted. The nuinbers associated with the lows embedded 
in the meanders and eddies along the coast correspond to 
temperature, not velocity. Note fully developed equatorward 
surface current with embedded meanders and eddies, along 
with poleward surface current north of y~1173 km. 
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North-South velocity at model day 426.0 


Figure 4.52 Experiment 5: Cross-section at y~704 km 
{41.25® N) of the meridional component of velocity at day 
426 (yd 60). Dashed lines indicate equatorward flow, while 
solid lines indicate poleward flow. Note strong poleward 
undercurrent along with sxibsurface flow associated with an 
eddy. 
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Flgura 4.53 Experiment 5: Surface velocity vectors at day 
426 (yd 60). To avoid clutter, velocity vectors are plotted 
at every third grid point in both the cross-shore and 
alongshore direction, and velocities less than 5 cm s'^ are 
not plotted. The numbers associated with the lows embedded 
in the meanders and eddies along the coast correspond to 
temperature, not velocity. 
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Temperature contour at model day 516.0 


Figure 4.54 Experiment 5: Surface ten 5 >erature at day 516 
(yd 150). The contour interval is 1® C. The temperature 
decreases towards the coast. 
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Temperature contour at model day 606.0 


Figure 4.55 Experiment 5: Surface tenperature at day 606 
(yd 240) . The contour interval if? 1® C. The temperature 
decreases towards the coast. Note maximum extent of 
upwelling. 
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Figure 4.58 Mid-monthly ECMWF winds for 1984 in m s*^; 

(a) January, (b) February, (c) March, (d) April, (e) May, 

(f) June, (g) July, (h) August, (i) September, (j) October, 

(k) November, and (1) December. The latitudinal grid point 
1 (129) corresponds to 35® N (47.5® N); so that the 
latitudinal grid point of 52 corresponds to 40® N. Maximum 
wind vector is 20 ro s~^. 
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Figure 4.58 (continued). 
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Figure 4,59 Experiment 6: Surface velocity vectors at day 
396 (yd 30). To avoid clutter, velocity vectors are plotted 
at every third grid point in both the cross-shore and 
alongshore direction, and velocities less than 5 cm s'^ are 
not plotted. The numbers associated with the low and high 
embedded in the meanders and eddies along the coast 
correspond to temperature, not velocity. Note equatorward 
surface flow with embedded meanders and eddies. 
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Distance off shore (km) 

-B53,3 -682.7 -512.• -341.3 -170.7 0.0 



North-South vel contour at model day 426.0 


Figure 4.60 Experiment 6: Surface meridional con^onent of 
velocity at day 426 (yd 60). The contour interval is 5 cm 
s* . The dashed lines Indicate eguatorward velocities. 
Latitudinal alongshore distance 0 (1408) km corresponds to 
35® N (47.5® N). Note weak poleward flow along the coast 
north of y-1273 km. 
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Figure 4.61 Experiment 6: Cross-section at 7-704 km 
(41.25® N) of the meridional component of velocity at day 
546 (yd 180). Dashed lines indicate eguatorward flow, while 
solid lines indicate poleward flow. Note existence of 
poleward undercurrent. 
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Temperature contour at model day 456.0 


Figure 4.62 Experiment 6: Surface tenperature at day 456 
(yd 90). The contour interval is 1® C. The temperature 
decreases towards the coast. Note start of upwelling. 
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Temperature contour at model day 666.0 


Figure 4.64 Experiment 6: Surface temperature at day 666 
(yd 300). The contour interval is 1® C. The temperature 
decreases towards the coast. Note occurrence of fllcunent at 
y~250 km. 
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Figure 4.66 Experiment 7: Total kinetic energy per unit 
mass time series for Climo III over the entire domain (Units 
of the kinetic energy are cm^s**) . Peak value of 166 cm^s*^ 
at day 915 (yd 185). 
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Figure 4.67 Experiment 7: Total kinetic energy per unit 
mass time series for Climo IV over the entire domain (Units 
of the kinetic energy are cm^s’^) . Peak value of 145 cm^s'^ 
at day 1293 (yd 198). 











Time (day) 


Figure 4.68 Experiment 7: Total kinetic energy per unit 
mass time series for Climo, Climo II, Cllmo III, and Climo 
IV over the entire domain (Units of the kinetic energy are 

cm^s'^) . 
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Figure 4.71 Experlirent 8: Total kinetic energy per unit 
mass time series for 1981, 1982, 1983, and 1984 over the 
entire domain (Units of the kinetic energy are cm^s"^) . 
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V. SnMMULT AND SSCGMMRNDATZONS 


A. SUMMARY 

The different experiment runs are summarized In Teddies 
4.2, 4.4, euid 4.5. In all first year mns (see Table 4.2), 
both a surface eguatorward current euid a poleward 
undercurrent developed. There was also a polemrard coastal 
surface current present during the winter season. Upwelllng 
was a seasonal feature which, as expected, had the coolest 
temperatures during a non-El Nlfio year (1981) emd warmest 
temperatures during an El Nlilo year (1983) . Both meanders 
and eddies were generated. The eddies were cold core, 
cyclonic, euad had a size around the Rossby radius of 
deformation (30 km). Along with an equatorward surface 
current, both mecuiders and eddies remained at the end of 
each year. Instability analyses revealed that both 
barotroplc euid bai-ocllnlc processes were importeuit eddy 
generation mechanisms. In all first year runs, no filaments 
were observed. 

The follow-on (second year) model runs are siumnarlzed In 
Ted)le 4.4. In general, a weeik coastal poleward surface 
current was evident during the winter season In the northern 
region of the domain for all years (Climo II, 1982 and 
1984). A mecuiderlng eguatorward surface current, already 
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present, at the start of these experiments was maintained. 

A poleward undercurrent was generated in all experiments. 

The eddies, already present at the start of the experiments, 
were both cold core, cyclonic, emd warm core, anticyclonic 
eddies in Climo II and 1982, while in 1984 there were only 
cold core, cyclonic eddies. The size and offshore extent of 
the eddies were consistent with the first experimental runs. 
Seasonal upwelling occurred with minimum ten^eratures of 
~12® C for Climo II emd 1982, euid of ~10® C for 1984. The 
year 1984 had the largest extent of cool upwelled water as 
expected for a year with the strongest equatorward winds. 
Filaments developed in all three e3q>eriments, but was most 
pronounced (as expected) in 1984. The filaments had cm 
offshore extent between 200 amd 350 km. Instadaility 
analyses revealed that both barotropic and baroclinic 
processes were involved in eddy generation. 

The last two years of the extended (four year) runs are 
summarized in Tadjle 4.5. Listed here are Climo IV, 1983, 
and 1984. Similar results are obtained: A poleward coastal 
surface current, a surface equatorward current, a poleward 
undercurrent, along with embedded meouidere amd eddies. Of 
interest, upwelling occurs later in 1983 (at yd 150) than 
for the other years (yd 60 to yd 90). Although the minimum 
temperature is the same (11.0® C) for all years, 1984 
experienced a larger area of cooler water (all along the 
coast) them that of the El Niflo (1983) year. In all 
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experiments, there were both cold and warm core eddies. 
Finally, filaments were observed during Climo IV and 19B4, 
but not during 1983. 

Experiments 1 through 6 investigated the formation or 
generation of features such as currents, meanders, eddies, 
and filaments. In Experiments 7 and 8 a (piasi-steady 
equilibrium state for the CCS was reached and the features 
observed to develop in previous years were maintained. 

Figure 5.1 plots the meain, eddy, and total kinetic 
energy per unit mass for Climo, 1981, and 1983 (Experiments 
1-3) together. During these first generation years. 

Figure 5.1a shows 1981 had the most total kinetic energy and 
1983 (El Nlho year) the least. This also holds true for 
both the mean emd eddy kinetic energy per unit mass plots 
(Figures 5.1b and 5.1c). The decrease in kinetic energy 
(total and mean) during the last few days of 1983 is not an 
error, but rather an adjustment for continuation with the 
start of the 1984 kinetic energy plots. 

Figure 5.2 plots the mean, eddy, cuid total kinetic 
energy per xinit mass for Climo II, 1982, emd 1984 
(Experiments 4-6) together (the second years of 
generation). The largest values of total, mean, and eddy 
kinetic energy were observed during 1984 (Figure 5.2a 
through 5.2c). Climo II had the least total and eddy 
kinetic energy (Figures 5.2a and 5.2c), but both Climo II 
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and 1982 had nearly identical low mean kinetic energy 
(Figure 5.2b). 

In all experiments, the atmospheric pressure cycle and 
subsequent movement of low and high pressure patterns 
affected the wind orientation and strength over the model 
domain. The surface current was observed to respond to the 
prevailing wind direction. As the wind turned eguatorward 
and exerted stress on the water surface, EkmcUi tremsport 
offshore resulted In the upwelling of cooler water along the 
coast. Interannual varleibllity was thus observed as 
expected. 

In conclusion, the experiments showed that, due to the 
wind forcing, features such as currents, meanders, eddies 
and filaments could be generated and maintained. These 
results support the hypothesis that wind forcing is a very 
In^ortant mechcuilsm for the generation of many of the 
observed features In the CCS. 

B. COMPARISON WITH OBSERVATIONS DURING AN EL NifiO YEAR 

It has been shown that the model results qualitatively 
agree with avalloUsle observations of the CCS during non-El 
Nliio years. For exaitple, the simulated surface and 
subsurface currents, along with eddies amd filaments agree 
In size and magnitude with available observations of the 
CCS, as listed In Batteen s!L (1989) . 
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Here the focus Is shifted to con^rlsons of observed and 
model results during an El Nlfio year. Qualitative results 
from both of the El Nlho experiments (3 and 8) along with 
some observations taken In the CCS during the 1983 El Nlfio 
are presented In Table 5.1. Because the model does not 
Include remote forcing, the oceanic teleconnectIon will not 
be expressed In the model results. There Is still good 
agreement between the model and observations tcdcen during 
the 1983 El Nlfio within the model domain. In particular, 
warmer sea surface teo^eratures along with onshore advectlon 
are evident, especially during the first three months of the 
year. Poleward surface flow (for example Figure 5.3) Is 
enhcuiced and the entire current system appears weaker theui 
In the other model years. Based on these results. It is 
concluded that atmospheric forcing played a exibstantlal role 
In the anomalous conditions observed In the CCS during the 
El Nlfio of 1983. 

C. NAVAL RELEVANCE 

As the Navy focus has recently shifted from blue water 
to the littoral, EBC models such as the one used here take 
on an Increased Isportance. Although this EBC model would 
require a finer scale before It could be used to predict the 
location of features of Interest, It Is a start in the right 
direction. 
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Accurately predicting the location of meanders, eddies 
and regions of upwelllng becomes critical for the successful 
use of mines, along with aviation, eutqphlblous, and covert 
swimmer operations. The meeuiders cuid eddies could Inqpact 
several aspects of mining (l.e., laying, known position, and 
mine removal). Smaller anphlblous craft could get swept off 
course If caught up In a strong jet associated with a 
meander. The cold upwelled water may be nutrient rich emd 
support a pleuikton bloom which could significantly reduce 
visibility In the water. In addition, the colder water 
adversely affects swimmer Immersion times. 

Fog often results In upwelllng regions where warmer air 
Is present over the cooler upwelled water. The effect of 
this fog Is to reduce surface visibility, which could be a 
plus (such as to provide cover for cm cutiphlblous Icmdlng) or 
minus (such as when defending cm area against covert 
operations). This decreased visibility may also affect 
aviation operations by requiring Instmment flight rules 
(IFR) vice visual flight rules (VFR). Finally, other 
optical emd Infrared devices will be degraded by the 
presence of Increased moisture In the atmosphere (l.e., 
fog) . 

D. SECGBSMEMDATZONS 

Future studies should let the model run for multiple 
years, to allow for spln-up. As seen In multiple year 
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kinetic energy per unit mass time series plots, it appears 
that spin-up may take only three to four years, but this 
needs to be confirmed. Additional studies should use daily 
rather than monthly winds to examine the effects of events 
and relcocations. The incorporation of bottom topography and 
a more realistic coastline would allow preferred eddy 
generation areas to be identified. The incorporation of 
bottom topography could also lead to an even better 
depiction of the undercurrent, which plays em importsmt role 
in eddy generation. Finally, finer horizontal amd vertical 
resolution should be incorporated into the model to aid in 
more accurately predicting the location of meanders, eddies, 
and fileunents. 


139 






TABLE 5.1 COMPARISON OF OBSERVED AND MODEL RESULTS DURING 
AN EL NINO YEAR (1983). 


El Nliio (1983) observations 
in CCS compared to meem 

Experiment 3 
(1983 cold 
restart) 

Eaqperiment 8 
(1983 extended 
run) 

Warmer SST 
(1. 2, 3) 

Simulated 

Simulated 

Wecdcer overall current 
system 
(2) 

Simulated 

Simulated 

Enhanced poleward flow 
(2, 3) 

Simulated 

Simulated 

Enhanced onshore advectlon 
(1, 2, 3) 

Simulated 

Simulated 

Maximum etjuatorward surface 
current farther offshore 
(2) 

Simulated 

Simulated 


References: (1) Rlenecker and Mooers (1986) 

(2) Slnqpson (1983) 

(3) Slnpson (1984) 
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Figure 5.1a Experiments 1-3: Total kinetic energy per unit 
mass time series for climatology (1980-1989), 1981, and 1983 
over the entire domain (Units of the kinetic energy are 
cm^s’^) . 
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Figure 5.1b Experiments 1-3: Mean kinetic energy per unit 
mass time series for climatology (1980-1989), 1981, and 1983 
over the entire domain (Units of the kinetic energy are 
cm^s'^) . 
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Figure 5.1c Experiments 1-3: Eddy kinetic energy per unit 
mass time series for climatology (1980-1989), 1981, and 1983 
over the entire domain (Unite of the kinetic energy are 
cm^s"^). 
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Figure 5.2a Experiments 4-6: Total kinetic energy per unit 
mass time series for Climo II, 1981, and 1983 over the 
entire domain (Units of the kinetic energy are cm^s"^) . 
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Figure 5.2b Experiments 4-6: Mean kinetic energy per unit 
mass time series for Climo II, 1981, and 1983 over the 
entire domain (Units of the kinetic energy are cm^s'^) . 
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Figure 5.2c Experiments 4-6: Eddy kinetic energy per unit 
mass time series for Climo II, 1981, and 1983 over the 
entire domain (Units of the kinetic energy are cm^s'^) . 
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North-South vel contour at model day 45.0 


Figure 5.3 Experiment 3: Surface meridional component of 
velocity at day 45. The contour interval is 1 cm s'^. Thi 
dashed lines indicate equatorward velocities. Latitudinal 
alongshore distance 0 (1408) km corresponds to 35® N (47.5 
N). Note enhanced poleward flow. 
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APPENDIX A. INSTABILITIES FOR EXPERIMENT 1 (CLIMO) 

This appendix contains additional cuialyses from 
Experiment 1 (climatological year) for the periods when the 
total kinetic energy per unit mass becemie quasi-steady. 
Energy transfer plots which consist of barotroplc (mesui 
kinetic energy to eddy kinetic energy) and barocllnlc (eddy 
potential energy to eddy kinetic energy) were produced for 
the following model days: 207 - 216, 228 - 231, cuad 297 - 
303. These plots appear as Figures A.l through A.6. 
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Figure A.l Experiment 1: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 
climatological (1980-1989) year, for model days 207 to 216 
The contour interval is 1 ergs cm*^ s'^. 
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Figure A.2 Experiment 1: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., baroclinic energy transfer) 
for climatological (1980-1989) year, for model days 207 to 
216. The contour interval is 1 ergs cm'^ s‘^. 
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Figure A.3 Experiment 1: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 
climatological (1980-1989) year, for model days 228 to 231. 
The contour interval is 1 ergs cm'^ s*^. 
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Figure A.4 Experiment 1: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., barocllnic energy transfer) 
for climatological (1980-1989) year, for model days 228 to 
231. The contour interval is l ergs cm'^ s'^. 
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Figure A.5 Experiment 1: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 
climatological (1980-1989) year, for model days 297 to 303. 
The contour interval is 1 ergs cm*^ s'^. 
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APPESIDIZ B. INSTABILITIES FOR EXPERIMENT 2 (1981) 


This appendix contains additional analyses from 
Experiment 2 (1981) for the periods when the total kinetic 
energy per unit mass became quasi-steady. Energy transfer 
plots which consist of barotropic (mean kinetic energy to 
eddy kinetic energy) and baroclinic (eddy potential energy 
to eddy kinetic energy) were produced for the following 
model days: 210 - 222, 240 - 252, 297 - 303, and 309 - 333. 
These plots appear as Figures B.l through B.8. 
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Figure B.l Experiment 2: Energy transfers of mean to eddy 
kinetic energy (i.e., barotroplc energy transfer) for 1981, 
for model days 210 to 222. The contour interval is 10 ergs 
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Figure B.2 Experiment 2; Energy transfers of eddy potential 
to eddy kinetic energy (i.e., baroclinic energy transfer) 
for 1981, for model days 210 to 222. The contour interval 
is 10 ergs cm'^ s'^. 
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Figure B.3 Experiment 2; Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1981, 
for model days 240 to 252. The contour interval is 10 ergs 
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Figure B.4 Experiment 2: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., baroclinic energy transfer) 
for 1981, for model days 240 to 252. The contour interval 
is 10 ergs cm'^ s'^. 
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Figure B.5 Experiment 2: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1981, 
for model days 297 to 303. The contour interval is 10 ergs 
cm’^ s'^. 
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Figure B.6 Experiment 2: Energy transfers of eddy potential 
to eddy kinetic energy (l.e., baroclinic energy transfer) 
for 1981, for model days 297 to 303. The contour Interval 
is 10 ergs cm'^ s*^. 
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Figure B.7 Experiment 2: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1981, 
for model days 309 to 333. The contour interval is 10 ergs 
cm'^ s'^. 
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Figure B.8 Experiment 2: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., baroclinic energy transfer) 
for 1981, for model days 309 to 333. The contour interval 
is 10 ergs cm'^ s‘^. 
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APPENDIX C. INSTABILITIES FOE EZPERIXENT 3 (1983) 


This appendix contains additional analyses from 
Experiment 3 (1983) for the periods when the total kinetic 
energy per unit mass became quasi-steady. Energy trcinsfer 
plots which consist of barotropic (mecui kinetic energy to 
eddy kinetic energy) and baroclinic (eddy potential energy 
to eddy kinetic energy) were produced for the following 
model days: 186 - 192, 222 - 234, 249 - 258, 279 - 288, 312 
- 315, 324 - 327, and 342 - 348. These plots appear as 
Figures C.l through C.14. 
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Figure C.l Experiment 3: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1983, 
for model days 186 to 192. The contour interval is 10 ergs 
cm‘^ s’^. 
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Figure C.2 Experiment 3: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., baroclinic energy transfer) 
for 1983, for model days 186 to 192. The contour interval 
is 10 ergs cm*^ s’^. 
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Figure C.3 Experiment 3: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1903, 
for model days 222 to 234. The contour Interval is 10 ergs 
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Figur* C.4 Experiment 3: Energy transfers of eddy potential 
to eddy Kinetic energy (i.e., baroclinic energy transfer) 
for 1983, for model days 222 to 234. The ccmtour interval 
is 10 ergs cm'^ s*^. 
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Figur* C.5 Experiment 3: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1983, 
for model days 249 to 258. The contour interval is 10 ergs 

F«flt '3 * • 1 


172 



Distance alongshore (km) 



Figure C.6 Experiment 3: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., baroclinic energy transfer) 
for 1983, for model days 249 to 258. The contour interval 
is 10 ergs cm*^ s'^. 
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Figure C.7 Experiment 3; Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1983, 
for model days 279 to 288. The contour interval is 10 ergs 
cm"^ s*^. 
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Flgur* C.8 Experiment 3: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., barocllnlc energy transfer) 
for 1983, for model days 279 to 288. The contour Interval 
is 10 ergs cm*^ s"^. 
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Figure C.9 Experiment 3: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1983, 

for model days 312 to 315. The contour interval is 1 ergs 
3 — * 1 
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Flgur* C.IO Experiment 3: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., barocllnic energy transfer) 
for 1983, for model days 312 to 315. The contour Interval 
Is 1 ergs cm'^ s'^. 
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Figure C.ll Experiment 3: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1983, 
for model days 324 to 327. The contour Interval is 1 ergs 
cm'^ s'^. 
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Figure C.12 Experiment 3: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., baroclinic energy transfer) 
for 1983, for model days 324 to 327. The contour interval 
is 1 ergs cm*^ s*^. 
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Figure C.13 Experiment 3: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1903, 
for model days 342 to 348. The contour interval is 1 ergs 
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Figure C.14 Experiment 3: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., baroclinic energy transfer) 
for 1983, for model days 342 to 348. The contour interval 
is 1 ergs cm*^ s‘^. 
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APPENDIX D. INSTABILITIES FOE EXPERIMENT 4 (CLIMO II) 


This appendix contains additional analyses from 
Experiment 4 (Climo II) for the periods when the total 
kinetic energy per unit mass became quasi-steady. Energy 
transfer plots which consist of barotropic (mecui kinetic 
energy to eddy kinetic energy) cuid baroclinic (eddy 
potential energy to eddy kinetic energy) were produced for 
the following model days: 378 - 408 (yd 12 - 42), 453 - 465 
(yd 87 - 99), 504 - 516 (yd 138 - 150), and 624 - 636 (yd 
258 - 270). These plots appear as Figures D.l through D.8. 
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Figure D.l Experiment 4: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for Climo 
II, for model days 378 to 408 (yd 12 - 42). The contour 
interval is 1 ergs cm'^ s"^. 
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Figure D.2 Experiment 4: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., baroclinic energy transfer) 
for Climo II, for model days 378 to 408 (yd 12 - 42). The 
contour interval is 1 ergs cm'^ s'^. 
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Figure D.3 Experiment 4; Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for Climo 
II, for model days 453 to 465 (yd 87 - 99). The contour 
interval is 1 ergs cm'^ s"^. 
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Plgur« D.4 Experiment 4: Energy transfers of eddy potential 
to eddy kinetic energy (l.e., barocllnlc energy transfer) 
for Cllmo II, for model days 453 to 465 (yd 87 - 99). The 
contour Interval Is i ergs cm*^ s'^. 
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Figure D.5 Experiment 4: Energy transfers of mean to eddy 
kinetic energy (l.e., barotropic energy transfer) for Climo 
II, for model days 504 to 516 (yd 138 « 150). The contour 
interval is 1 ergs cm'^ s*^. 
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Figure D.6 Experiment 4: Energy transfers of eddy potential 
to eddy kinetic energy (l.e., baroclinic energy transfer) 
for Cllmo II, for model days 504 to 516 (yd 130 - 150). The 
contour interval is 1 ergs cm'^ s*^. 
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Figure D,7 Experiment 4: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for Cllmo 
II, for model days 624 to 636 (yd 258 - 270) . The contour 
interval is 1 ergs cm'^ s"^. 
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Figure 0.8 Experiment 4: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., barocllnic energy transfer) 
for Climo II, for model days 624 to 636 (yd 258 - 270). The 
contour Interval is 1 ergs cm'^ s*^. 
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APPEIIDZZ B. ZNSTABZLITZBS FOB BZPERZKBHT 5 (1982) 


This appendix contains additional cuialyses from 
Esiperiment 5 (1982) for the periods when the total kinetic 
energy per unit mass became quasi-steady. Energy treuasfer 
plots which consist of barotropic (mecin kinetic energy to 
eddy kinetic energy) euad baroclinic (eddy potential energy 
to eddy kinetic energy) were produced for the following 
model days; 507 - 516 (yd 141 - 150), and 525 - 528 (yd 159 
- 162). These plots appear as Figures E.l through E.4. 
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Figure B.l Experiment 5: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1982, 
for model days 507 to 516 (yd 141 - 150). The contour 
interval is 1 ergs cm'^ s‘^. 
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Figure E.2 Experiment 5: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., baroclinic energy transfer) 
for 1982, for model days 507 to 516 (yd 141 - 150). The 
contour interval is 1 ergs cm'^ s'^. 
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Figure B.3 Experiment 5: Energy transfers of mean to eddy 
kinetic energy (l.e., barotropic energy transfer) for 1982, 
for model days 525 to 528 (yd 159 - 162). The contour 
Interval is 1 ergs cm'^ s'^. 
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Figure B.4 Experiment 5: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., baroclinic energy transfer) 
for 1982, for model days 525 to 528 (yd 159 - 162). The 
contour Interval is 1 ergs cm’^ s'^. 
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APPENDIX F. ZNSTABZLZTZBS FOR BZPBRZMBNT 6 (1984) 

This appendix contains additional emalyses from 
Experiment 6 (1984) for the periods when the total kinetic 
energy per unit mass became quasi-steady. Energy transfer 
plots which consist of barotroplc (mecui kinetic energy to 
eddy kinetic energy) and barocllnlc (eddy potential energy 
tc eddy kinetic energy) were produced for the following 
model days: 576 - 591 (yd 210 - 225), 606 - 624 (yd 240 - 
258), emd 645 - 657 (yd 279 - 291). These plots appear as 
Figures F.l through F.6. 
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Figure P.l Experiment 6: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1984, 
for model days 576 to 591 ^d 210 - 225). The contour 
interval is 10 ergs cm'^ s" . 
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Figure F.2 Experiment 6: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., barocllnic energy transfer) 
for 1984, for model days 576 to 591 (yd 210 - 225). The 
contour interval is 10 ergs cm"^ s’^. 
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Figure F.3 Experiment 6: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1904, 
for model days 606 to 624 (yd 240 - 258). The contour 
interval is 10 ergs cm"^ s"^. 
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Figure F.4 Experiment 6: Energy transfers of eddy potential 
to eddy kinetic energy (i.e., barocllnic energy transfer) 
for 1984, for model days 606 to 624 (yd 240 - 258) . The 
contour interval is 10 ergs cm"^ s'^. 
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Figure P.5 Experiment 6: Energy transfers of mean to eddy 
kinetic energy (i.e., barotropic energy transfer) for 1984, 
for model days 645 to 657 (yd 279 - 291) . The contour 
interval is 10 ergs cm'^ s" . 
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